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ABSTRACT 


A  criterion  is  developed  for  bubble  initiation  from  a  gas  filled 
cavity  on  a  surface  in  contact  with  a  superheated  layer  of  liquid.  It 
is  found  that  the  temperature  of  bubble  initiation  on  a  given  surface  is 
a  function  of  the  temperature  conditions  in  the  liquid  surrounding  the 
cavity  as  well  as  the  surface  properties  themselves.  It  is  also  found 
that  the  delay  time  between  bubbles  is  a  function  of  the  bulk  liquid 
temperature  and  the  wall  superheat  and  is  not  constant  for  a  given 
surface. 

By  consideration  of  the  transient  conduction  into  a  layer  of 
liquid  on  the  surface,  a  thermal  layer  thickness  is  obtained.  With 
this  thickness  and  a  critical  wall  superheat  relation  for  the  cavity, 
a  bubble  growth  rate  is  obtained. 

Bubble  departure  is  considered  and  it  is  found  that  the  Jakob 
and  Fritz  relation  works  as  long  as  the  true  (non-equilibrium)  bubble 
contact  angle  is  used.  The  effect  on  the  departure  size  of  the  virtual 
mass  in  the  surrounding  liquid  is  found  to  be  negligible  at  one  gravity. 
That  is,  at  one  gravity  the  primary  effect  of  bubble  growth  velocity 
on  bubble  departure  size  is  found  to  be  due  to  contact  angle  changes. 

The  initiation,  growth  and  departure  criterions  are  each  experi¬ 
mentally,  individually,  checked.  They  are  then  used  to  compute  the 
heat  transfer  near  the  knee  of  the  boiling  curve  using  only  an  experi-  • 
mental  determination  of  the  number  of  bubbles  as  a  function  of  wall 
superheat  and  other  known  quantities.  Finally  the  q  vs.  Tw  -  Tgat 
relation  is  computed  and  measured  and  compared.  The  comparison  is 
satisfactory. 
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Introduction 

When  the  wall  temperature  exceeds  the  saturation  tempera¬ 
ture  of  the  fluid  which  is  in  contact  with,  a  thin  layer  of  super¬ 
heated  fluid  near  the  wall  is  formed.  If  there  is  a  cavity  on  the 
wall,  initially  filled  with  inert  gas,  a  bubble  will  start  to  grow 
from  that  spot,  when  the  wall  superheat  becomes  sufficiently 
high.  A  further  increase  in  the  wall  superheat  will  cause  an 
increase  in  the  growth  rate  and  bubble  generation  frequency. 

The  growth  of  the  bubble  will  lift  up  the  superheated  liquid 
layer  from  heating  surface.  The  departure  of  the  bubble  will 
carry  away  the  thermal  layer  from  an  influence  circle  around 
the  nucleate  site.  This  repeated  process  gives  rise  to  the  high 
heat  transfer  in  nucleate  boiling. 

When  the  bubble  frequency  exceeds  a  certain  value,  the 
distance  between  consecutive  bubbles  is  so  small  that  they  join 
together  into  an  unstable  and  shaky  chain  and  the  idealizations 
made  in  the  analytical  portions  of  this  study  will  no  longer  be 
valid  for  an  accurate  prediction  of  the  heat  transfer  rate. 
Reference  (7)  gives  the  boundary  dividing  the  isolated  bubble 
region  from  the  columns  of  bubbles  region.  This  work,  in 
general  applies  only  to  the  region  of  isolated  bubbles. 
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1. 


Bubble  Initiation  Theory 
a.  General  Description 


A  bubble  is  generally  initiated  from  a  small  gas  filled  cavity 
or  crack  on  a  solid  surface  so  long  as  the  surrounding  fluid 
is  heated  to  a  sufficiently  high  temperature.  Both  the  pre- 
existence  of  the  gas  phase  and  the  temperature  are  necessary 
but  not  sufficient.  This  mechanism  has  been  completely  dis¬ 
cussed  in  Reference  (3),  in  the  case  of  homogeneous  tempera¬ 
ture  field.  An  extension  of  this  mechanism  to  the  non- 
homogeneous  temperature  field  (which  is  of  primary  interest) 
will  be  developed  in  this  section.  A  similar  situation  can  be 
seen  in  the  process  of  chicken  incubation.  A  chicken  can  be 
incubated  from  an  egg,  if  and  only  if,  this  egg  is  a  fertile  one. 
The  bubble  initiation  mechanism  is  similar  in  that,  a  gas  phase 
must  already  exist  along  with  the  right  temperature  conditions 
as  shown  in  Figure  1. 
b.  Transient  Thermal  Layer  . 

Since  the  convection  intensity  near  a  solid  wall  is  damped  down 
due  to  the  no  slip  boundary  condition  for  a  solid  surface,  the 
use  of  the  pure  conduction  equation  is  justified  in  determining 
the  temperature  distribution  in  this  thin  layer  of  fluid  near  the 
heating  surface.  For  this  particular  problem,  a  simplified 
physical  model  is  shown  in  Fig.  2. 

Initial  condition  is 


T  =  Tw  x—°  -j 

T  =  T„  *+  x  vo  J 

4 

Boundary  condition  is 

T=s  Tw  tf-f  | 

T  =  Too  X  *=■  °o  J 

The  solution  to  this  problem  is 


±  =  o 

~t  ^  o 

found  from  Ref.  (1)  as 


(1) 

(2) 


T~T„  =  (Tw-T„>  *-fc  ^ 

X* 

‘d'T  _  _  Tw ~ Tqg 

2*  spritt 


(3) 

(4) 


2 


at  x  =■  o 

2L 

>eo 


(- 


.  _  TW-Tcp 

'  y/irkt 


(5) 


(6) 


If  the  actual  temperature  distribution  near  the  wall  is 
assumed  to  be  a  straight  line  distribution,  the  slope  of 
this  straight  line  is  determined  by  equation  (5).  This 
assumption  has  been  justified  through  measurements  made 
in  Reference  (6).  With  this  fact,  one  can  introduce  the 
notion  of  thickness  of  transient  thermal  layer  by  drawing 
a  tangent  line  from  x  =  0  on  the  T  -  T^^  X  curve  defined 
by  (3),  the  interception  of  this  straight  line  on  X-axis  gives 
the  transient  thermal  layer  thickness  as  shown  in  Fig.  (3). 

S  =  Jnr  1l± 

This  means  that  the  temperature  distribution  at  any  instant 
varies  linearly  from  the  wall  to  x  =  $  ,  beyond  £  ,  the  fluid 
does  not  know  whether  the  wall  is  hot  or  cold.  The  layer 
thickness  increases  with  the  square  root  of  waiting  time, 
c.  Criterion  of  Bubble  Growth  Initiation 

Having  the  defintion  of  the  transient  thermal  layer,  one  can 
determine  the  time  required  from  the  beginning  of  generation 
of  thermal  layer  to  the  beginning  of  bubble  growth.  This 
period  is  defined  as  the  waiting  period  of  a  bubble,  t  .  The 
criterion  for  initiating  bubble  growth  is  from  Reference  (4) 
that  the  thermal  layer  surrounding  the  bubble  nucleus  must  be 
at  a  mean  temperature  equal  to  or  above  the  temperature  of 
the  vapor  in  the  bubble  in  order  to  give  rise  to  an  inward  flow 
of  heat  from  the  superheated  thermal  layer  to  the  bubble 
through  the  bubble  wall.  Before  bubble  growth  the  bubble  is 
in  the  condition  of  thermo-static  equilibrium.  The  equation 
of  static  equilibrium  for  the  bubble  is  then 

Ap  =  zg- 

(7) 

With  the  help  of  the  Clausius  -  Clapeyron  thermodynamic 
equilibrium  relation,  one  has 

L 


Ap  = 


At 


t  JL 

Uf  fir  f 


ATjjL 

Tut 


(8) 
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where 


^  T  =  ~Tb  “  Ts„f 

^  P  Jb  ' 

~Tb  are  temperature  and  pressure  of  the  vapor 

in  the  bubble  at  the  initial  stage  of  bubble 
growth. 

Eliminating  from  above  equations  yields 

2.^  7i#+ 


^  T 


"Ti-  ~  ^Sat 


or 


~Tb~  7~oo  ~  71  »  7^0 


2  p-  Tsm¥ 


(9) 


During  the  waiting  period,  the  bubble  wall  can  be  treated 
approximately  as  an  insulated  hemi-spherical  surface  of 
radius  Rc  (the  temperature  distribution  in  the  surface 
tension  layer  of  the  bubble  is  unknown).  Presumably  there 
is  tangential  conduction  in  a  thin  layer  around  the  bubble  so 
that  the  interface  temperature  is  constant.  A  physical  model 
of  this  idealization  is  shown  in  Fig.  4. 


From  potential  flow  theory  and  the  fluid  flow  analogy,  the 
potential  line  in  fluid  flow  is  just  equivalent  to  the  isothermal 
line  in  heat  conduction,  the  distance  of  an  isothermal  line 

3 

passing  through  the  top  point  of  a  waiting  bubble  is  Rq 
distant  from  heating  surface  when  measured  on  the  straight 
part  of  this  isothermal  line. 


Fluid  temperature  at  X  =  — 2“Rc  *s 

^  =  (Ttv-rm)(i-ij^)+-ro0=  Tw-trw-Tm) 


BRc 
2  S 


(10) 


Equating  this  temperature  to  the  bubble  temperature  yields 
the  criterion  of  initiation  of  a  bubble  growth  from  a  nucleate 
site  of  cavity  radius  Rc  as 


or 


c  Ju 

£ _  3  (~7jy  ~~  T«o)  Pc 


3KC 

ZS 


Jw  ~  Tsmt  (  1  Rcf  u  ) 


(ID 
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when  £  is  expressed  in  terms  of  the  waiting  period 


wft  +mk 


r  ( r„  ~~  To* )  He 

L  7w~  TSq+  ( 1  + 


z—\ 

ref0u) 


The  Most  Favorable  Cavity  Radius  for  Initiating  Bubble 
Growth  and  the  Minimum  Waiting  Period _ 

As  the  waiting  time  increases,  the  thermal  layer  increaes 
until  to  a  certain  condition  such  that  the  temperature  line 
of  luid  becomes  tangent  to  the  bubble  temperature  curve. 
At  this  instant,  if  and  only  if,  there  is  a  dry  cavity  of 
radius  R  f  on  the  heating  surface,  a  bubble  will  begin  to 
grow  from  this  spot.  This  radius  corresponding  to  a 
minimum  waiting  period,  is  called  the  most  favorable 
cavity  radius  Let  us  now  turn  our  attention  to  the 

solution  of  equation  (11). 

Solving  for  Rc  from  (11)  yields 


S(TW-T,.^  r  ±  ] 

3  (Tw-Tmi  1  »  (Tw-T**fSPL  J 


For  any  given  waiting  period,  there  are  two  possible  cavity 
radii  which  will  nucleate.  When  these  two  cavity  radii  are 
equal,  it  means  the  two  intersecting  points  coincide,  (see 
Fig.  5),  or  the  fluid  temperature  line  and  bubble  tempera¬ 
ture  curve  are  tangent  to  each  other.  Observing  (13)  gives 
the  condition  of  equal  roots  of  Rc  as 

I  _  /Z  (  Tty  ~  Tt p )  Ism-t  d  ^ 

(Tw-Ts<+)*SfiL  ~ 

Solving  for  $  which  is  by  definition  $min  yields 


/z  (~7yy~~  Tap)  Ts*+ 
ft  L  tTw~TSm+)z 


Equal  root  condition  in  (13)  with  help  of  (14)  gives 

KCr  =  ^_)  __  4-Ts*+  & 

T  3  (Tw-llo)  )(1 

( -t„)  =  =,  j±±t&znofJkxt. d*  (fcji 

'TM.  ryrk  fp'L*(7w-TM)  + 

Upper  and  Lower  Bounds  of  Radius  of  Active  Nucleate  Cavib 
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The  thermal  layer  cannot,  in  general,  increase  without 
limit  with  the  waiting  time.  It  will  be  washed  off  by  the 
natural  convection  of  the  fluid  as  it  grows  beyond  the  thick¬ 
ness  of  natural  convection  layer  he  This  means 
f  _  C 

°m«x  °nc  /j 

Knowing  Smmx  from  the  natural  convection  information, 
the  maximum  and  minimum  cavity  radius  for  initiating  a 
bubble  growth  can  be  calculated  from  (13)  with  help  of  (16) 


W4 


max i  min. 


r,  +  zEjgl  1 


(17) 


f. 


Any  cavity  outside  this  interval  cannot  qualify  as  an  active 
nucleation  site.  A  diagram  is  shown  in  Fig.  5. 

A  Numerical  Example  for  a  Quantitative  Illustration 
Fluid  =  degased,  distilled  water 
TSit+  =  212°  F  =  672°  R 
Too  =  202°  F 
Tv*  =  242°  F 
O'  =38.3  10"4  lb.  /ft. 
fu  =  0.  0374  lbm/  ft.3 
'k  =  1.82  10"6  ft.  2/sec. 

L  =  755.  10 3  ft.  lb.  /  lbm. 


The  most  favorable  cavity  radius  for  initiating  a  bubble  growth 
is  from  (IS)  a 


■4  &  Tsm+ 


/.  2/ 


/•*  ft 


7  f„L(Tw-Ts„) 

The  corresponding  thermal  layer  thickness  is  from  (14) 

C  _  Izo'  Tsg+  (  Tw  -  7k2  —  4  o  t  .  ,2^  n. 

*">■"  ~  f.  L  (Tw-TMft  *  97  "  ft 

The  minimum  waiting  period  for  this  bubble  is 

S'* 

j.  \  _  Omin 

wJm;»  nr  ft 


=  4 .  /  -4* 


•  /  o  *  Sec  . 
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2. 


Bubble  Growth  Theory 
a*  Assumptions 


i.  Neglect  any  convection,  not  due  to  bubble  itself, 
completely. 

ii.  Neglect  the  change  of  mass  of  fluid  due  to  evaporation 

or  condensation  through  the  interface  of  vapor  and  liquid. 

iii.  One  dimensional  case  is  converted  into  the  three  dimen¬ 
sional  case  by  the  introduction  of  a  curvature  factor. 

iv.  Neglect  the  inertia  force  and  the  surface  tension  of  the 
fluid. 

v.  Constant  properties  of  fluid. 

vi,  Spherical  bubble  surface. 

vii.  Uniform  wall  temperature,  mainfluid  temperature  and 
fluid  pressure. 

After  a  waiting  period  t  ,  the  bubble  is  going  to  grow. 

For  the  first  few  moments,  the  surface  tension  effects  and 
the  inertia  effects  of  surrounding  fluid  are  so  large  such 
that  the  growth  rate  is  controlled  by  momentum  equation 
but  after  the  radius  increases  to  about  twice  its  initial  value, 
the  surface  tension  and  inertia  effects  will  become  negligible, 
so  that  the  growth  rate  is  controlled  only  by  the  heat  transfer. 
In  this  study,  only  the  heat  transfer  effects  will  be  considered 
for  the  evaluation  of  bubble  growth  curve, 
b.  Formulation  and  Solution 

For  simplicity,  the  one  dimensional  physical  model  for  heat 
transfer  mechanism  is  given  in  Fig.  6. 

Initial  condition  is 


T  =  Tkv  -*TW-Tm )  ^  for 

where  §  _ 


T—  To 


for 


Boundary  condition  is 

T =  Tsa+  j-or 

T=  Toe  -h' 

Introducing  a  new  variable 


o  XL*  g 


l  x  -c.  ao 


J 


~t  ZAO 


*  =°  }  *  >« 

X  a*  J 


(18) 


(19) 
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(20) 


0  ==■  T"  Tu+  such  that 

0^  —  —  Tlv^"7saf  v  ^oo  ~  Too-Tso* 


then  (18)  and  (19)  are  transformed  to 

A 


0 : 

1 

« 

t 

n 

~fw~Too  - 

f 

& 

i 

11 

TSo+ 

=  &«, 

S' 

a 

=  o 

■for  x 

& 

II 

8~* 

1 

■TSat 

c  ^or  o<x*$ 

for  5C  ^  § 


y 


a  o 

(21) 


1 

■for  ^  =  os  j 


-t  >o 


(22) 


The  governing  equation  is  then 

M  _ _L  2L 

&  9t  <23) 

The  problem  is  now  reduced  to  a  semi-infinite  conductor, 
with  a  prescribed  initial  temperature  B(x,0)  =  and 

surface  temperature  zero,  then  the  solution  of  (23)  with 

conditions  (21)  and  (22)  will  be,  from  Reference  (1) 

j*>  r  tx-x'f  Cx+x'J* 

I  I  r.  .....  f  ~  4  k-t  ”  ***  ) 

-  e  j 


6  = 


2jnrkt 


r  ]%’  f 


(24) 


where 


99 

ax 


=  Q - fyy  ^  x  -for  o  <  pc  S 

u”  $  J 

for  f  -c  x'^oo^ 
r"w  (xtx'j 

[(.*.*>  dx  (25) 


“f  CX  y  -  |/^y 

=  0, 


at  x  =■  o  (25)  becomes 
99  \ _ / 


f— ) 

V  4jnr{M 


—  *'* 

x" ;  (2  x^)  e 


X ' 

~$£k  .  > 

e 


/ 


(26) 
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Referring  to  the  bubble  growth  model  as  shown  in  Fig.  8, 
the  governing  equation  for  bubble  growth  is 


dR 

or  df 


%  ft  A±r(2l)  +  % 

?•  al 


(27) 


where  ft>^  =  curvature  factor  where  I 

9 9  =  surface  factor  =  21TK  __  i^cosg 


-  base  factor  =  SHE*  Ujisl.  = 


Sin 


J- 


(28) 


_  _  = 

t>  -^nrf^y  <4' 

ft^  =  volume  factor  —  ~  Tl27r*3Q-CoSfj]’t’  3  n*  s/nfc+sj 

2  +  Cqs  f(  2+S/hp) 


4-rrr  R* 
3 


ft  =  contact  angle 

Jl  -  coefficient  of  heat  transfer  from  heating 

surface  to  the  steam  bubble  through  its  base  area. 
Substituting  (26)  into  (27)  yields 

=  ?<?s  4tcf  (f)  _  CTTW_J_)  +  3 

%  ftL  -r-****  ywJ  %fuL 

. (29) 

For  the  case  of  a  bubble  growing  in  an  infinite  fluid  field  of 
superheat  9W  then  %  -1-  %  - x*  %-*•  00  * 

and  (29)  becomes 

dp_  =  ?cAc/>  gu  _g_  ip 

d*  AL  JtFZF  pr  fuL  ft  (30) 

From  homogenions  solution  of  bubble  growth  rate,  one  has 
from  Scriven's  theory  (see  Reference  (18)) 

jn_  =  [t  e„fc  nr 

df  «/V~  ftL  J-t  (31) 

Comparing  (30)  and  (31),  one  finds  the  value  of  curvature 
factor  to  be 


%  for  J 


S  R 


(32) 


9 


Another  extreme  Case  is  for  J  -  nr,  it  reduces  exactly  to 
one  dimensional  case  then 

fc  ~ 1  (33) 

For  -  o  ,  and  it  reduces  to  Plesset’s  thin  layer 

case  which  gives 

Tl  —  (34) 

Combining  these  three  extreme  cases,  one  can  manufacture 
a  such  that  it  satisfies  (32),  (33),  and  (34)  simultaneously, 

/  r\  O  /oc\ 


where  Ff  is  the  time  average  of  bubble  radius  or 

ft  =  J-  f'Kdt  (36) 

Integrating  (29),  with  respect  to  time  t  gives 

/V 

t?_p  —  %  fc  f  I  [/)  -ID  -6  1  erJ-I—  IdJ-  + 

K  *e  nz J.jr&  P*  **  **  *  1  T»t>L 

or 

r,  r  jr  II  t  tflL  er/J= 

*  *  wmrJ  *  ^ 1  ^ 

*  #  ^  - 1  "Z1  jAt  ']  *  -ffi”  *  <3,) 


Normalizing  (37)  by  introducing  dimensionless  variables 
T  =  -"--t  ^  Jv  =  -y-  leads  to 


l-l  =3Z'fre»  fjE-  *#/?«' 

%  f,L  L’l'jr  J 

-  2  ~rfc  %l±Sa>-  -c 

J  ^  J  */&#.*•  (38) 


A  bubble  growth  plot  for  >*7^  ,  7^,  —  7^^  and 

y  y  ^  in  normalized  coordinates  is  shown  in  Fig.  (10). 
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c. 


Experimental  Result  Compared  with  the  Theory 


Fluid  =  Distilled  and  degased  water 


C 

V 

r 

f 


1.807  10'6  ft.2/ sec. 


1.007  btu/lbm. 

3.16  lO-6  ft.  2/sec. 

10-4  //p. 

38.3  lO-4  lb. /ft. 

59.97  lbm/ft.  3 
fv  =  0.0374  lbm/  ft.  3 
L  =  755  103  ft.  lb.  /lbm 
P  =  1  atm 

Surface  =  16  k  gold  polished  on  clothe  wheel  by  No.  8 
diamond  compound 
J*  -  0.  750  radian 

Data  recorded  - 
=  229.98°  F 
Toe  =  205.02°  F 
212°  F 

Bubble  number  1  = 


Camera  speed  =  1140  frames  /  sec. 
“tw  =  0.  0245  sec. 

=  0.  0166  sec. 

=  0.397  •  10~2  ft. 

§  =  0.372  •  10‘3  ft.  (from  (6)) 

K=  0.01097  •  10"3  ft.  (  from  ( 13)) 

For  =  1.  52 

fu 


5* 10"3  sec. 

R  =  3.  36  • 

o 

i 

03 

10  *  10~3  sec. 

3.89  • 

10"3  ft. 

15  •  10“^  sec. 

4.03  * 

10~3  ft. 

From  (37) 


with  A  =  0 

u 
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Bubble  number  2  = 

Camera  speed  =  1260  -frame  *  /sec. 

~tw  =  0.  0437  sec. 
f4  =  0.  0167  sec. 

Rd  =  0.  533  •  10"2  ft. 

$  -  0.498  *  10-3  ft.  (from  (6)  ) 


For 


%% 

% 


1.62 


=  5  * 

CO 

1 

o 

rH 

sec. 

R  =  3.99  * 

10‘3  ft. 

10  • 

10'3 

sec. 

4.93  • 

10"3  ft. 

15  • 

10'3 

sec. 

5.  25  • 

10'3  ft. 

From  (37) 
with  ^  =  0 


Bubble  number  3  = 

Camera  speed  =  1380  frame  /sec. 
=  0.  0275  sec. 

-  0.  0145  sec. 

Rd=  0.47f-  /O3 -ft. 

$  -  ■  io*  ft. 

For  =  1.73 

% 


=  5  * 

10“3  sec. 

R  =  3.928  * 

10“3  ft. 

10  • 

10“3  sec. 

4. 619  * 

10"3  ft. 

15  • 

10“^  sec. 

4.848  * 

10"3  ft. 

From  (37) 
with  ^  =  0 


The  corresponding  bubble  growth  curves  with  a  comparison  from 
theoretical  ones  are  shown  in  Fig.  11  The  bubble  growth  history 
for  these  three  bubbles  are  shown  in  Fig.  12,  13,  and  14. 


d.  Discussion  of  the  Bubble  Growth  Theory 

In  the  bubble  growth  theory,  the  thermal  layer  on  the  bubble 
surface  is  assumed  to  be  picked  up  by  a  growth  of  bubble 
immediately  at  the  last  moment  of  waiting  period.  From  the 
high  speed  photographic  study  described  above,  one  can  see 
that  at  the  first  moments,  the  bubble  growth  rate  is  very  high 
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and  the  bubble  expands  laterally  at  such  a  rate  that  in  fact 
a  very  large  piece  of  thermal  layer  is  picked  up  during  the 
first  few  moments.  This  fact  will  give  a  strong  support  of 
the  one  dimensional  approach.  Actually  the  bubble  growth 
history  is  composed  of  three  periods,  namely  the  waiting 
period  “tw  ,  the  unbinding  period  ~tub  ,  and  the  departure 
period  f  .  When  the  wall  superheat  increases,  the  waiting 
period  of  bubble  at  a  particular  cavity  decreases  very  rapidly. 
If  the  thermal  layer  thickness  calculation  is  still  based  on 
the  waiting  period,  the  error  will  be  very  large.  This  will 
make  the  deviation  between  the  theoretical  bubble  growth  rate 
and  actual  one,  very  large.  From  Fig.  15,  in  which  the 
dynamic  effect  and  surface  tension  effects  to  bubble  growth 
are  shown,  the  following  can  be  seen:  During  the  waiting 
period,  the  bubble  is  heated  in  order  to  initiate  growth  from 
its  cavity.  During  unbinding  period,  the  bubble  is  trying  to 
librate  itself  from  the  binding  force  of  surface  tension  and  the 
inertia  effects  of  its  surrounding  fluid.  The  bubble  radius 
increases  very  slowly  and  the  momentum  equation  governs  the 
motion  of  bubble  surface.  During  the  departure  period,  the 
effects  of  surface  tension  and  inertia  of  fluid  become  so  small 
that  the  heat  transfer  equation  governs  the  motion  of  bubble 
surface  and  the  thermal  layer  is  picked  up  by  the  growing 
bubble  immediately  during  the  first  few  moments  of  this 
period.  Therefore  the  thermal  layer  thickness  for  very  high 
wall  super  heat  case  where  the  waiting  period  is  very  short 
should  be  calculated  by 

%  —  J^7rA  instead  of  %  ~  . 

Observations  from  those  22  bubbles  listed  in  subsection  2-c 
show  that  the  departure  period  was  nearly  constant,  the 
waiting  period  changed  by  a  factor  eight  to  one.  The  uneven 
heating  due  to  a  500  watt  light  source  for  photography  purpose 
at  the  rear  side  of  test  section  caused  a  pronounced  unsymmet- 
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rical  turbulent  convection  of  fluid  which  changed  the 
thermal  layer  distribution.  The  temperature  fluctuations 
associated  with  the  turbulence  gave  rise  to  fluctuations  in 
the  waiting  period  also . 

Departure  Criterion 
a.  Formulation 

From  the  bubble  growth  equation  and  the  normalized  bubble 
growth  diagram  in  Fig.  10,  one  can  see  that  a  bubble  can 
either  depart  from  its  nucleate  site,  stay  there  or  collapse 
there.  For  the  case  of  Too  ^  ,  the  ubble  grows  mono- 

tonically,  the  bubble  must  eventually  depart  from  the  heating 
surface  due  to  monotonically  increasing  buoyant  force  of  the 
bubble.  For  the  case  of  Too  <  TSot  ,  it  is  not  clear,  one 
needs  a  criterion  for  judging  if  or  not  a  bubble  will  depart 
from  the  heating  surface.  In  case  where  the  bubble  departs, 
the  time  to  departure  is  the  quantity  of  interest.  This  is 
the  next  question  to  be  considered. 

In  order  to  study  the  dynamical  departure  criterion  of  a 
bubble  from  a  heating  surface,  a  force  and  motion  analysis 
is  necessary.  The  assumption  of  a  perfectly  spherical 
bubble  will  give  no  information.  For  this  reason,  one  needs 
to  modify  the  physical  model  of  bubble  from  a  spherical  one 
to  some  other  shape.  In  this  section,  the  inertia  force  effect 
of  the  surrounding  fluid  is  calculated  by  truncated  spherical 
bubble  model  which  is  of  course  not  exact.  After  writing 
down  the  governing  equations,  having  chosen  the  dominating 
variables  of  the  quasi -static  solution  which  is  compatible  . 
with  the  analytical  solution  in  the  static  case,  the  assumption 
of  a  spherical  bubble  for  evaluating  bubble  volume  and  surface 
is  resumed. 

From  potential  flow  theory,  the  inertia  mass  of  the  surrounding 
fluid  of  a  solid  sphere  departing  from  a  solid  plane  boundary 
is  (Reference  2 ) 


(39) 


By  the  fact  of  non-sphericity  of  the  bubble  shape  and  with 
support  of  experiments,  Davidson  had  corrected  this  value 

t0  =  -a- ft? 

where  iP  is  the  volume  of  bubble,  f  is  the  density  of  the 
fluid. 

M*  is  assumed  to  be  uniformly  distributed  as  a  very  thin 
mass  layer  on  the  surface  of  the  bubble.  Considering  the 
bubble  as  a  thin  shell  loaded  with  hydrostatic  force  and  iner¬ 
tia  force,  defining  p  as  the  pressure  on  the  inner  face  of 

'  TP 

surface  tension  layer  of  the  bubble  and  as  the  pressure 

on  the  outer  face  surface  tension  layer  of  the  bubble,  and  with 
the  concept  of  inertia  mass  layer,  one  has  by  momentum 
equation  from  Fig.  16. 


z.,:  riving) 


R.J-  s..,- ?](}.-}>+ rt 


(40) 


(41) 


k'}  ~  fcti}  (42) 

Shell  formula  of  force  equilibrium  applied  at  the  top  point  of 
the  bubble  gives 


at  point  (  X  ,  j  ),  it  becomes 


S  fr)<p 


Defining  the  equation  of  meridian  curve  of  the  bubble  as 


J  J  C  X,  b  s  -t) 

where  b  is  a  function  of  time,  then 

/+(&)*_  S;„ 


p- 

ctx 


K  = 


dx * 


(43) 


(44) 


(45) 


(46) 


Substituting  (40),  (41),  and  (42)  into  (43)  and  (44)  and  elimi- 
«_  /  / 

nating  Y  s  and  p  s  yields 


15 


Zrf"  C 1  Sin<p  O'  fp  _p  \  -  ht  ) 

~b  *  ~~T~  ')  I*'* 

—  <w;j+4{*  ^  <47> 

Using  an  approximation 

if  =  S  =  4-nrK  (48) 

in  which  R  is  determined  by  (37)  and  is  function  of  time 
substituting  (46)  into  (47),  the  differential  equation  for  the 
bubble  surface  can  be  found. 

The  initial  conditions  are  that 

Q-t  -t=o,  b~Rc  ,  }  =  Hc  (49) 

A  numerical  method  similar  to  the  technical  of  Bashforth 


and  Adams  in  Reference  13  is  necessary  for  a  complete 
solution. 


where  -  angle  of  contact 

Therefore  from  solution  of  (47),  putting  f-%  ,  one  has 

b,  -t>  (51) 

where  aT \  and  £  are  functions  of  time 

O 

or  *■.  0.,  b,  ir)  (52) 

Then  the  departure  criterion  is 


X~°  (53) 

from  which  the  departure  time  can  be  solved  by 

*•  t}!v '  -t]=  o  (54) 

If  (54)  has  no  real  positive  root  t^,  it  means  no  departure. 


If  (54)  has  a  real  positive  root  t^  andj(  t^)^o  f  it  means 
that  the  bubble  will  depart  at  t^.  This  is  a  complete  descrip¬ 


tion  of  formulation  and  method  for  finding  out  an  exact  solution. 
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b.  Solution 


An  approximate  solution  for  this  problem  can  be  obtained 
in  the  following  way. 

For  case  of  very  small  contact  angle  %'<<^  >  (50)  can 
be  simplified  by  assuming  a  nearly  spherical  bubble  and 
Xa  smf9  negligible  compared  with  ^  in  the  last  term 
of  (50). 

[’ - <*«*>]  (55) 

The  law  of  motion  with  help  of  the  notion  of  the  inertia  mass 
layer  gives 

fP  2prx° «  SMf. 

With  help  of  (44)  and  (41),  the  above  equation  is  reduced  to 

fP fy-  wi ~  (~&  *  | > +*prx.e si»fm 

Sin  f<,  _  J51  ) 


or 


'/-//"/  -  >  ■*  (■ 


Adding  a  term  —  C  k  +  R R)  on  both  sides  of  above 


equation  gives 


+  W  -  Zik  ( & 


*  (  ^Sin(Po _ 51 


Ro 


) 


For  small  contact  angle  %-  l  /^,  2 
one  has  then 


A- nr 


*  , 


v-r^n<  -  is  -  Vf  ^  -  fj* 

. ' . (56) 


Now  dropping  out  all  subscripts  0  in  (55)  and  (56)  and 
keeping  their  original  definitions  yields 

X  +  Uks.  =JL-!mu  r  _ ue+Rk)} 

K  X  b  o'  i-  4S<f-fu\ JR  J 


nf 


(r 


1 


=  mx  S'*?.  -$r)+  "■'£&  Ur+RR) 
^  js?  -4-8 


(57) 
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For  the  static  case  -Tt  so  ,  all  dot  and  double  dot  terms 
in  (57)  drop  out,  so  (57)  is  reduced  to 

/  <ss;nf  z 

x  ~  b  tr 


A/ 

H 

V 


(a) 


m*  o' 


f  S:r»q>  _L  \ 
</-/»>/  v  *  “  » 


(b) 


} 


(58) 


The  departure  criterion  for  the  governing  equation  (58)  was 
given  in  Reference  (9)  by  Fritz.  He  used  the  numerical 
analysis  result  of  Bashforth  and  Adams  (Reference  13)  to 
correlate  the  dimensionless  quantity  — against  the 

the  contact  angle.  He  concluded  that  the*bubble  will  depart 
at  maximum  volume,  so  that 


A  nrxzd  /  SV/t  f  _  J_  } 

**  departure 


3  3 

=  Q3I3Q 
s 


? 


(59) 


where 


2IT 


(so) 

which  is  called  Laplacian  characterizing  length  of  a  static 
bubble  (59)  is  supported  successfully  by  experiments  for 
f  =0  up  to  -  2. 4  radians . 

Comparing  (57)  and  (58)  term  by  term,  one  can  see  that  a 
quasi-static  solution  exists,  if  the  static  Laplacian  length 
is  replaced  by  a  dynamic  length  such  that 


a 


=  a.  [/-  — 11  ^ - (kl+*ik) I  2 

4  *-  J 


48i. t-pjk  j  (6D 

then  the  quasi-static  criterion  for  dynamic  departure  will  be 

I?  [o.  313  a  f+  jq&'u  K+  **>%-] 

or4 


3  ~3 


t par4ure. 


0.313  Qg  J 


I- 


nf 

24R 


,  .*  ...  aj  / 

(-4  R  +R  HJ  departure 


(62) 


where  ^  is  the  dynamic  contact  angle. 

This  departure  criterion  is  based  on  an  assumption  which  gives 
the  right  result  in  the  static  case  and  yields  an  estimate  of  the 
effect  of  dynamic  forces  on  the  departure  size  of  that  bubble. 
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Crudely  one  can  regard  the  dynamic  forces  as  altering  the 
gravity  that  the  bubble  sees.  Observations  from  the  high 
speed  photographs  of  bubbles  show  that  the  contact  angle 
apparently  changes  with  the  velocity  of  the  triple  interface. 

It  is  experimentally  well  established  that  the  contact  angle 
only  has  its  equilibrium  value  when  the  system  is  static. 

In  our  experiments  (see  Fig.  12,  13  and  14)  the  change  in 
the  contact  angle  was  obvious  as  the  bubble  went  from 
advancing  to  receding  and  also  when  the  rate  of  growth 
changed.  The  departure  size  was  found  to  be  a  function  of  the 
receding  dynamic  contact  angle  rather  than  the  mean  value 
experienced  by  the  bubble  which  is  in  contact  with  the  surface. 
The  most  important  force  change,  the  contact  angle,  appears 
to  be  the  viscous  force.  To  a  first  approximation,  the  change 
in  contact  angle  might  be  linear  in  the  ratio  of  the  viscosity 
force  and  surface  tension  force.  Dividing  bubble  Weber  num¬ 
ber  with  bubble  Reynolds  number  yields  a  new  number. 


which  is  the  ratio  of  viscosity  force  to  the  surface  tension 
force.  (See  the  appendix) 

The  correlation  relating  J  to  is  then 

J  =  (64) 

where  A  is  a  constant  to  be  determined  by  experiments. 

In  Reference  (10),  Staniszewski  had  performed  51  experi¬ 
ments  using  water  and  alcohol  as  the  fluids  at  different 
system  pressure,  A  can  be  evaluated  using  these  experi¬ 
ments  and  is 

A  =  6ZSO  (65) 

Actual  measurements  of  the  dynamic  contact  angle  from 

Figs.  12,  13  and  14,  show  that  (65)  was  in  good  agreement. 

Since  Staniszewski1  s  experiments  were  done  at  relatively  low 
wall  superheat,  the  dynamical  effects  other  than  the  contact 
angle  change  are  secondary  to  the  gravity  effects  on  departure. 
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This,  we  feel,  is  the  primary  cause  of  apparent  dynamic 
effect  on  departure  size  by  Staniszewski.  The  value  of 
contact  angle  reported  by  Staniszewski  are  actually  aver¬ 
ages  for  a  large  number  of  readings.  It  is  assumed  here, 
on  the  basis  of  our  own  observations,  that  the  dependence  of 
the  departure  size  on  the  bubble  growth  rate  is  a  result  of 
a  change  in  contact  angle  due  to  dynamic  effects  rather  than 
any  dynamic  effects.  In  the  subsequent  subsection,  a  cal¬ 
culation  will  be  made  to  show  that  the  dynamic  effects  in  the 
liquid  are,  indeed  negligible  for  the  condition  under  which 
our  data  and  that  of  Staniszewski' s  were  taken  . 

(62)  is  only  an  approximate  criterion  for  the  case  of  small 
contact  angle.  The  exact  one  should  be  found  from  (53) 
which  is  too  complicated  to  work  out  for  this  study  as  far  as 
the  time  is  concerned.  However  with  the  suggestion  re¬ 
ported  in  subsection  3  -  a.  ,  and  a  complete  understanding  of 
Reference  13,  a  general  departure  criterion  might  be  worked 
out  with  the  help  of  machine  computation. 

Putting  J  —  0  in  (47),  a  criterion  for  zero  gravity  case  can 
be  obtained. 

To  show  the  dynamic  contact  angle  effects  on  bubble  departure, 
Stanizewski's  data  is  shown  in  Fig.  18  in  which 

=  0.4*.*  | 

*<>  =  J  (66) 

where  R.  is  in  ft. 

is  in  radian 
is  in  ft.  /sec. 

c.  The  Period  to  Departure  and  the  Bubble  Generation  Frequency 
Putting  =^R  andR=R(j,  one  can  solve  for  R^  from  (62). 
Having  solved  for  R^,  the  corresponding  time  t^,  the  departure 
period,  can  be  solved  by  (37). 
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d. 


With  t ,  found  from  (37)  and  t  found  from  (12)  the  bubble 
d  w 

generating  frequency  is  then 


f  =.  _ L _ 

A  bubble  generating  cycle  diagram  is  shown  in  Fig. 


Discussion 


In  this  analysis,  the  effect  of  the  disturbance  of  the  surrounding 
fluid  due  to  natural  convection  is  completely  discarded.  Ac¬ 
tually  the  natural  convection  of  the  surrounding  fluid,  the 
irregularity  of  bubble  shape,  the  surface  condition  of  the  wall, 
the  disturbances  arising  from  a  growth  and  a  departure  of  the 
neighboring  bubbles,  the  bubble  population  density  will  strongly 
influence  the  departure  diameter.  Therefore,  a  deviation  of 
only  +  10%  of  the  departure  diameter  from  the  experimental 
result  is  not  surprising 
Comparison  with  Experimental  Results 
Fluid  =  Distilled  and  degased  water 

Surface  =  No.  8  diamond  compound  polished  gold  surface 
f  =  0.  750  radian 

Data  are  exactly  the  same  as  that  in  bubble  growth  theory. 


For  bubble  number  1  - 

*R,  =  4.009  *  10  ^  ft.  from  (62)  and  (48) 

d  -3 

=  3.974  •  10  ft.  from  experiment 

For  bubble  number  2  - 

R,  =  5.363  ‘  10'3  ft.  from  (62)  and  (48) 

d  -3 

R^  =  5.328  *  10  ft.  from  experiment 

For  bubble  number  3  - 

R,  =  4.886  •  10~3  ft.  from  (62)  and  (48) 

d  -3 

R^  =  4.  792  *  10  ft.  from  experiment 

*  The  value  R  in  <p  was  taken  from  the  slope  of  experimen¬ 
tal  bubble  growth  curve  at  t  =  t^  insiead  of  iheore+ical  Ones. 
A  plot  of  measured  departure  radii  against  the  calculated  ones 
is  shown  in  Fig.  19. 
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4. 


Heat  Transfer  Correlation 


a.  Explanation  of  Boiling  Curve 

Boiling  curve  can  be  best  explained  by  the  theory  of  "bulk 
convection  of  the  transient  thermal  layer".  Observations 
show  that  when  the  wall  temperature  exceeds  the  saturation 
temperature  of  the  fluid,  the  heat  transfer  increases  very 
rapidly  with  the  wall  temperature.  Many  researchers  have 
tried  to  explain  why  this  occurs.  The  following  study  explains 
these  observations  by  means  of  a  so-called  theory  of  bulk 
convection  of  the  transient  thermallayer,  or  simply  bulk  con¬ 
vection  theory.  When  the  boiling  starts,  the  bubbles  depart 
from  the  heating  surface.  In  departing,  the  bubbles  bring 
part  of  the  layer  of  superheated  liquid  adjoining  the  bubble 
into  the  main  body  of  fluid.  At  the  same  time,  the  cold  fluid 
flows  onto  the  heating  surface.  The  heat  transfer  rate  for  the 
first  few  moments  after  this  process  is  very  high  due  to  the 
very  high  temperature  gradient  near  the  wall.  After  a  certain 
time,  a  new  thermal  layer  is  built-up,  and  a  new  bubble  starts 
to  grow.  When  this  bubble  grows  to  a  certain  size,  it  departs 
from  the  heating  surface  and  a  new  thermal  layer  is  brought 
to  the  main  body  of  fluid  again.  By  this  kind  of  repeated  trans¬ 
portation  of  thermal  layer  (which  is  technically  called  bulk 
convection),  heat  is  transferred  to  the  fluid  from  the  wall.  The 
heat  transfer  rate  by  this  process  is  nearly  proportional  to  the 
square  root  of  bubble  generation  frequency.  In  Fig.  20,  one 
can  see  that  from  A  to  B,  heat  transfer  rate  increases  very 
rapidly  due  to  the  increase  in  T  -  T  ,  which  increases  the 
bubble  generating  frequency,  the  enthalpy  content  of  the  tran¬ 
sient  thermal  layer  and  the  density  of  active  cavity  population. 
At  B  the  active  cavity  population  has  been  increased  to  a  sat¬ 
uration  state  such  that  the  influence  circle  of  each  bubble 

touches  one  another.  A  further  increase  of  T  -  T  ,  does 

w  sat 

not  increase  area  of  production  of  transient  thermal  layer,  but 
the  bubble  frequency  and  enthalpy  content  of  thermal  layer 
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continues  to  increase.  Therefore  after  B  the  rate  of  increase 
of  J  *s  reduced.  B  is  a  point  of  inflection.  From  B  to  C 
the  bubble  frequency  increases  until  to  a  certain  stage  such 
that  unstable  and  shaky  vapor  jets  are  formed.  These  con¬ 
tinuous  vapor  columns  reduce  the  effective  area  of  production 
of  transient  thermal  layer,  such  that  the  curve  becomes  con¬ 
cave  downward.  From  C  to  D,  the  effective  area  of  production 
of  transient  thermal  layer  decreases  more  rapidly  than  in¬ 
crease  of  the  enthalpy  content  in  the  thermal  layer  due  to 
increase  of  T  -  T  ,  therefore  the  curve  drops.  At  point 

W  S  al 

D,  the  effective  area  of  production  of  transient  thermal  layer 
has  been  reduced  to  zero,  a  steady  and  continuous  blanket  of 
vapor  exists  between  the  heating  surface  and  main  fluid.  The 
fluid  gets  essentially  no  chance  to  touch  the  heating  surface, 
therefore  no  transient  thermal  can  be  built  up  on  the  heating 
surface  and  the  heat  transfer  rate  reaches  to  a  minimum 
value.  Bulk  convection  process  is  completely  stopped  at  D. 

A  further  increase  of  T  -  T  ,  will  increase  heat  flux  again 

W  S  al 

by  radiation  and  conduction  across  the  gap. 
b.  Mechanism  of  Heat  Transfer 

The  heating  surface  in  pool  boiling  is  divided  into  two  parts, 
the  bulk  convection  area  and  the  natural  convection  area.  In 
the  area  of  bulk  convection,  heat  is  assumed  to  be  transferred 
into  the  fluid  by  transient  conduction  process.  Following  the 
departure  of  a  bubble  from  the  heating  surface,  a  piece  super¬ 
heated  liquid  is  brought  into  the  main  body  of  the  fluid.  By 
this  kind  of  repeated  process  heat  is  transferred  from  heating 
surface  to  the  main  body  of  the  fluid.  In  the  area  of  the  natural 
convection,  heat  is  supposed  to  be  transferred  from  heating 
surface  into  the  main  body  of  fluid  by  the  usual  convection  pro¬ 
cess  in  a  continuous  manner.  A  physical  model  of  bulk  con¬ 
vection  mechanism  is  shown  in  Fig.  21  . 
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At  stage  1,  a  piece  of  superheated  transient  thermal  layer 
is  torn  off  from  heating  surface  by  the  departing  bubble, 
and  at  the  same  time,  the  cold  fluid  from  the  main  body  of 
the  fluid  flows  onto  the  heating  surface,  after  a  time  interval 
t  ,  this  cold  liquid  layer  is  heated  to  a  condition  such  that 
the  tiny  bubble  on  that  cavity  is  able  to  grow  which  is  shown 
as  stage  2.  At  stage  3,  the  bubble  grows  laterally  with  a 
very  high  rate  such  that  a  very  large  piece  of  thermal  layer 
is  picked  up  in  a  very  short  time  interval.  At  stage  4,  the 
bubble  is  going  to  depart  from  the  heating  surface  which  will 
bring  the  situation  immediately  to  stage  1  again.  This  cyclic 
process  furnishes  a  way  to  transfer  the  heat  from  the  heating 
surface  to  the  main  body  of  the  fluid. 


c. 


The  system  which  is  used  to  evaluate  the  heat  transfer  per 
bubble  cycle  is  as  follows. 

Formulation 

i,  Natural  Convection  Component 

A  theoretical  study  of  natural  pool  convection  points  out 
that  the  natural  convection  heat  transfer  can  be  corre¬ 


lated  by  using  two  dimensionless  groups  namely 

'V' 

The  Nusselt  Number  ^  O 

No  =  -r 

fck. 

The  Rayleigh  Number  =  /;  (tw-t„)d 

■h  p 


For  Laminar  Range 

/  O**  2  •  /  0  i  ] 

A/u  =  O.S4-  R.Z.  [ 

For  Turbulent  Range 

2  I07  ^  3- 10  ^ 

M/  =  o.  /4  Ra  / 


where  D  = 
A  = 


Area  of  heating  surface. 


(68) 


(69) 
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This  correlation  was  first  studied  experimentally  by 
cryder  and  Finalborgo  and  was  summarized  by  Fishenden 
and  Saunders. 


Substituting  (68)  into  (69)  and  making  use  of  the  definition 
of  heat  transfer  coefficient  yield 
For  Laminar  Range 


(70) 


/o*  <  R*  <  2  10 

kc-£frw-Ti>- 

For  Turbulent  Range  f 

Z  IO*  R*  3  !0  4, a  .3 

U  -  It  (  Tvy-Tc.)  =  o-'+fc  J  (VI) 

For  illustrative  purpose,  a  numerical  example  is  given 

as  follows. 

Liquid  =  water 
,o 


/  =  59.97  lbm/ft.  3 

=  242”  F  C  =  1. 007  btu/lbm  °F 

o 


T  .  =  212"  F 
sat 

rp  _  n  /i  nO 

w 


=  202"  F  3  =32.2  ft.  / sec. 

=  1  7/8"  =  0.  156  ff 
=  1.81  •  10"6  ft.  2/sec. 


T* 

D 

r 

>  =  0.316  I O* 


=  10"4  V°F 

-tr 


From  (68)  ^  ?!  =  $t;4  lo  >  2.  ■  !o? 

his 

which  is  in  turbulent  range,  from  (71) 

]yc  =  o  /^c  [22l2kz2k2l£lf=  1.7  B+u/(ft  SecJ 

The  thickness  of  the  thermal  layer  of  natural  convection  is 


S. 


fck. 


(72) 


Wc  —  (Tw-Too)  =  o.2t+-/S*fl 

ii.  Bulk  Convection  Component 

From  equation  (2),  one  can  obtain  the  heat  transferred 
through  unit  area  of  heating  surface  to  the  fluid  during 
time  t  as 

jrr-Tjcj>dx=cj><T„-T^  £Z£SfM(73) 
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For  this  case.  S'  is  not  a  constant  throughout  the 
bubble  base  where  the  transient  conduction  thermal 
layer  is  developing.  Such  a  donut -shaped  layer  is 
illustrated  in  Fig.  22. 


For  convenience  in  integration,  the  initial  state  is 
taken  at  the  end  of  waiting  period,  so  that 

S  ^ 

i  =  J/Trb+J  m  Sw  ( 

Making  use  of  (72),  the  heat  transferred  into  transient 
thermal  layer  as  well  as  in  the  main  body  of  fluid  beyond 
the  transient  layer  during  one  bubble  formation  cycle  is 

aO=  f  *- 1 2 -nr dr)  +  nrO?;-ti) 

J*‘  Hi 

=  2pc(Tw-T<~>  [  [  27 rrS*l-  -h  nr  (&*- R  x)  £  ] 

or  ^  y  ^ 


where  is  influence  radius 

R,  =  2  R,  for  the  isolated  bubble  case  ^ 

1  d  \  (76) 

R  <  2R^  for  the  close  packed  case  J 

Since  Rc  <<  R^,  and  §*  is  nearly  linear  in  T  ,  so 
(75)  can  be  approximated  to  yield 

=  Zfc  (  tw-Tc )  f  *  * Sd- i  fC  (  5*  - (77) 

If  "n"  is  the  number  of  active  cavities  of  radius  Rc  per 
unit  area  of  heating  surface,  and  -j-  is  the  frequency  of 
bubble  generation,  then  the  heat  transfer  rate  per  unit 
area  due  to  bulk  convection  of  the  transient  thermal  layer 
is  approximately  from  (77) 

?.c—/  *<?  -  *fc<Tw-Tm>»f  F^](78) 

iii.  Vapor  Convection  Component 

In  addition  to  the  heat  transferred  directly  to  fluid,  heat 
is  also  transferred  directly  into  bubble  through  the  heating 
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surface  exposed  to  the  vapor  inside  the  bubble.  This 
component  is  important  only  at  very  high  wall  super¬ 
heat.  It  can  be  omitted  in  nucleate  boiling  region  but 
must  be  considered  in  the  film  boiling  region. 

Sf»? 

?,C  =  M (Tw-T,'t)h,  +  J 


(79) 


where  =  R  S7n^>  =  base  circle  radius  of  bubble 

Ac  =  Surface  area  of  cavity 

^  =  Heat  transfer  coefficient  of  vapor  con¬ 

vection. 

Since  Ac  and  Rc  are  very  small  quantities,  (79)  can  be 
reduced  roughly  to 

/?; s/4 


J  '  --  “  -  J  (80) 

iv.  General  Expression  of  Heat  Transfer 

Combining  (70)  (or  (71)),  (79)  and  (80)  leads  to 

} =  bc  *  b c  +  /•*- 

=  ( /-  nrn  £*)  (Tw-Too)  + 

-+  X  77 ^  ~  (81) 

A  three  dimensional  sketch  of  the  heat  transfer  as  a 
function  of  subcooling  and  wall  superheat  is  given  in 
Fig.  23  in  which  the  effects  of  subcooling  of  main  fluid  and 
wall  superheat  can  be  easily  interpreted  by  means  of  the 
bulk  convection  theory. 

Discussion 

The  population  density  of  bubbles  at  the  close  packed  condition 
is  such  that  the  bubbles  are  so  densely  packed  that  the  influence 
circle  of  one  nucleate  cell  touches  its  neighbors,  considering 
one  half  cell  as  indicated  in  Fig.  24  by  shaded  area,  one  has 

A/c  p  _  .2~  _  / _ _ 


72 cp  ~  A  ~  i  (2ft.J3H.-J  ~  <82> 

where  R^  =  2  R^  (83) 

(83)  was  justified  by  some  rough  experiments  in  which  a  ball 
a  was  pulled  up  from  the  bottom  of  water  tank 


of  radius 


which  has  a  layer  of  chalk  powder  on  the  bottom.  Observations 
showed  that  the  chalk  powder  within  a  circle  of  radius  R.  »  2  a 
moved  toward  the  center  forming  a  vortex  ring  in  the  wake  part 
of  the  ball.  This  vortex  ring  is  a  method  of  scavenging  away 
the  thermal  layer  within  this  influence  circle  and  putting  a  new 
layer  of  cold  liquid  on  the  heating  surface  bounded  by  the  in¬ 
fluence  circle.  A  sketch  of  this  process  is  shown  in  Fig.  25. 


From  potential  flow  theory,  the  velocity  potential  and  stream 
line  function  in  the  surrounding  fluid  of  a  departing  sphere  of 
radius  a  from  a  solid  plane  boundary  are  from  Reference  (2) 


#  =  -J£  ( 

2  \  [3 


+  a*_L+  MlC- 


a 


4h* 

(D  =  -  2JL  (-  *3*  +  ctll  +  a  V 
v  -M  t  Ah* 


Sh'r1 


_  _  a  ~  _  a  Z 

3*h‘  +h*r  3zh*r 


Tti*  -)c>sB  +  Cl  ( 8 •*) 

•  )  siie+c'  (9$) 


The  velocity  components  at  point  P  in  radial  and  meridian 
directions  are 


_  1 I L 13 

~  Tl  /•» 


+  £+£  -Q*  -JL-  - 
4f,3  32 h6  32 ft* 


)Cosfr 


(86) 


10  _ 

"O  ”  t  M  "  *  1 


6 


.  SZ  ,  RL-  +£ _ o 

f*  '  4h3  3*h*  Shsr*  Hk€t*  '  ,nd  (87) 


h  =  in  (86)  and  (8  7)  gives  roughly  the 


putting  a  =  R 

velocity  components  of  the  fluid  surrounding  a  departing 
bubble  as  r  r  ,  jp* 

Vr  *  +-  ' 

**  "i?  +  ^  *!n$  +" 

where  U  is  the  bubble  rising  velocity  at  departure  . 


(88) 


Example  with  experimental  data  taken  from  Reference  (5). 

^  =  8.3  •  10^  btu/ft.  2 /hr.  =  23  btu/ft.  ^/sec. 

T  =  230°F  evaluated  from  bubble  growth  theory  based  on 
the  bubble  growth  curve  given  in  Reference  (5). 

Tsat=  212°  F 
Tw  =  196°  F 

t  =  12  *  10  ^  sec. 

w 
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t  ■,  =  12  *  10  ^  sec. 

a 

R,  =  0.040  in  =  0.0033  ft. 
a 


2.6^. 


From  (74) _ 

f  =  T^TFrf  =  */•*  cyc/e/sec. 

Fig.  5  in  Reference  5  shows  that  this  is  the  close  packed 


10*jt 

3.70 


case 


R^  =  2  R^  =  0.0066  ft.,  (82)  gives 

K  *=*  4 —  „  =  £6  3  o 


From  (78),  the  heat  transfer  due  to  bulk  convection  is 


which  is  about  75%  of  the  total  heat  flux  supplied  to  the  heating 
surface. 


The  heat  flux  required  for  evaporation  of  vapor  into  bubble  is 

»/ L  it**/  -  i-* 1  B*a/f  &ecJ 

where  /5  _  itSZT  t? 3 

u»>ax  3  n4 

which  is  about  6%  of  the  total  heat  flux. 


The  heat  transfer  due  to  natural  convection  for  close  packed 
case  is  from  (71)  and  (82)  . 

>  0./4.P"  & 

■  J3  J 


L  ■  0- 


P 


2  J- 

-r- 


0./3  2 


rj*-smv 


which  is  about  0.  6%  of  the  total  heat  flux. 

Most  of  the  difference  between  the  calculated  and  experimental 
results  is  due  to  geometric  idealization.  From  above  example, 
one  can  see  that  the  bulk  convection  of  transient  thermal  layer 
from  heating  surface  to  the  main  fluid  constitutes  a  chief  means 
of  heat  transfer.  Bubbling  is  the  only  natural  mechanical  driving 
force  which  propels  such  a  bulk  convection.  The  bubble  growth 
theory  and  the  departure  criterion  in  nucleate  boiling  heat  trans¬ 
fer  are  important  not  because  they  can  carry  a  large  amount  of 
heat  due  to  evaporation  of  fluid  into  the  bubble  voids  (only  a  few 
percent^  but  because  they  supply  the  way  to  take  off  the  transient 
thermal  layer  repeatedly  from  the  hqating  surface. 
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e.  Comparison  with  Experiment  Results 

Two  sets  of  experimental  results  are  presented  in  this  section; 
one  set  is  our$^,  while  the  other  was  taken  from  Ref.  17, 
both  results  were  compared  with  this  theory. 

i,  Result  1  -  (From  our  experiments) 

Fluid  used:  Distilled  degased  water 

Surface:  Gold  layer  plated  on  copper  base,polished 
with  No.  8  diamond  compound. 

System  pressure  P  =  | 

Data  point  1  .  =  0.0620  btu/sec.  from  (91) 

T  =  218.73°  F 

w 

Tgat=  212. 00°  F 
T«>  =  178.56°  F 

N  =  12. 

N.  =  12  of  R  =  3.0460  *  10"5  ft. 
i  c 

from  (17),  (^c)mjLn  was  taken  as 

the  cavity  radius,  since  (R  ) 

J  c  max 

is  nearly  a  hundred  times  larger 

than  the  surface  texture  dimension 

N  =  0 
a 

Qp  =  0.0620  btu/sec.  from  (81) 

Data  point  2  =  Qp  =  0.  1202  btu/sec.  from  (91) 

T  =  235.09°  F 
w 

Tgat=  212.00°  F 
Too  =  199.72°  F 
N  =  18. 

N  =  12  of  R  =  3.  0460  •  10“5  ft. 
a  c 

f  =  69.  15  /sec.  from  (17). 

(11),  (12),  (62), 
(37),  and  (67) 

N,  =  6  of  R  =  0.7859  •  10-5  ft. 
i  c 

R  ,  =  4.15-  10'3  ft.  from  (62), 

(37) 

Qp  =  0.  1142  btu/sec.  from  (81) 
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Data  point  3  = 


Data  point  4  = 


Data  point  5  = 


Qp  =  0.  1433  btu  /sec. 

T  =  237. 11°  F 
w 

T  =  212°  F 
sat 

T,,,  =  201.87°  F 
N  =  20. 

IN  =  18.  .12  of  R  =  3.046  •  10-5  ft. 

^  •  I  Q 

i  f  =  78.46  1/sec. 

.  '6  of  R  =  0.  7859  •  10~5  ft. 

c 

f  =  53.08  1/sec. 

N  =  2  of  R  =  0.  7240  *  10"5  ft. 
i  c 


Rj  =  4.215  '  10  3  ft. 
d 

Qp  =  0.  1412  btu  /sec. 


Qp  =  0.  1866  btu/sec. 

T  =  237. 61°  F 
w 

T  =  212°  F 
sat 

T„  =  201.38°  F 
N  =20 


(N_  =  20.  ,12  of  Rc  =  3.046  •  10"5  ft. 

f  =  80.72  1 /sec. 


6  of  R  =  0.  7859  *  10“5  ft 
c 

f  =  61. 56  l/sec. 

„  2  of  R  =  0.  7240  *  IQ'5  ft 
c 

f  =  6. 44  1  /sec. 


*N,  =  0 

R,  =  4.231  '  10"3  ft. 
d 

Qp  =  0.  1584  btu/sec. 


Qp  =  0.  2157  btu/sec. 

T  =240. 65°  F 
w 

T  =212.  00°  F 
sat 

Tm  =  200.53  F 
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N  =  20. 


fNa 


=  20. 


< 


12  of  Rc  =  3.046  •  10"5  ft. 

f  =88.03  1/sec. 

-  6  of  R  =  0.  7859  *  10”5  ft. 
c 

f  =  87.06  1/sec. 

^  2  of  R  =  0.7240-  10~5  ft. 
f  =  78 . 60  1/sec. 


IN  =  0 

R  ,  =  4.322  *  10"3  ft. 
d 

Qp  =  0.2056  btu  /sec. 

A  comparison  of  experimental  result  with  theoretical  ones  is 
shown  in  Fig.  26. 

ii,  Result  2  -  (From  Fig.  8  on  Reference  (17)  ) 

Fluid:  n  -  pentane  C,-  H-  9 


Surface:  Nickel,  4/o  polished 

System  pressure:  1  atm 

Properties  of  fluid: 

T  =  97°  F 
sat  „ 

f  =37.8  lbm/ft.  a 
fu  =  0.  187  lbm/ft.  3 
ff"  =9.79-  10"4  lb/ft. 

^  =  146  btu/lbm  =  1136  •  10_3ft.  lb/lbm 

P  =  4.41  •  10'6  ft.  2 /sec. 

' k '  =  1.097  •  10~6  ft.  2/sec. 

c  =  0.527  btu/lbm 
)  =  8.1  •  10"3  /p*  (From  Ref.  (17)) 

Data  point  1  =  |R  .  390  btu  /(ft>  2  gec_ } 

T  =  111°  F 
w 

Tsaf  T«  =  97  F 
n  =  =  430  1/ft.2 

ina  =  ° 

\nt  =  430  of  Rc  =  0.  365  *  10“5  ft.  from(17) 


|p  =  0.390  btu  /(ft.  2  sec. ) 
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Data  point  2  = 


Data  point  3  = 


Data  point  4  = 


h 

T 

w 

Tsat 

n 


=  0.865  btu/(ft.2  sec.) 

=  119°  F 
=  =  97°  F 

=  1580  1/ft.  2 

=  430  1/ft.  2  of  Rc  =  0.  365  *  10_5ft. 

f  =  39.8  1/sec. 

from  (17), (11),  (12), (62),  (37), 
and  (67) 

=  1150  1/ft.  2  of  Rc  =  0.  233  •  10"5ft. 


from  (17) 

Rd  =2.60*  10"3ft.  from  (62)and(37) 
=  0. 806  btu  /  (ft.  2  sec. ) 


?R 

T 

w 

Tsat 


=  1.  208  btu 
=  122°  F 


/(ft.  2 
97°F 


sec. ) 


n  =  2480; 

(n  =  1580.  430  of  R  =  0.  365  •  10~5  ft. 

<  f  =  43.0  1/sec. 

<  1150  of  Rc  =  0.  233  *  10"5ft. 

f  =  42.  6  1/sec. 

n.  =  900  of  R  =  0.  205  *  10"5  ft. 
l  c 

Rd  =  2.  79  *  10~3  ft. 

^  =  1258  btu  /(ft.  2  sec. ) 


1500  btu  /(ft.  2  sec. ) 
124°  F 

=  97°F 


n  =  3800: 
(na  =  2480 


4 


430  of  R  =  0.365  *  10'5ft. 
c 

f  =  45. 4  1  /sec. 

1150  of  R  =  0.  233  *  10“5ft. 
.  c 

f  =  45.2  1/sec. 
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Data  point  5  = 


= 


R 


l900  of  R  =  0.  205  *  10~5ft, 
c 

f  =  44.  9  1/sec. 

1320  of  R  =  0.  1903  *  10~6  ft. 
c 

283  ‘  10"3  ft. 


^.p  =1611  btu  /(ft.  3  sec. ) 


% 


R 


Tw  - 


n 


1. 815  btu  /(ft.  Z  sec. ) 
5760. 


n  = 


v-5 


3800  ,  430  of  R  =  0.  365  *  10  "ft. 
c 


f  =  47.4  1/sec. 


1150  of  R  =  0.  233  ’  10 
c 


-5 


ft. 


f  =  47.3  1/sec. 


.-5 


900  of  R  =  0.  205  *  10  "ft. 
c 


f  =  47.  1  1/sec. 


v-5. 


n,  =  1960 


1320  of  R  =  0.  1903  *  10  "ft. 
c 


f  =  46. 3  l/sec.  ) 


-3 


R  ,  =  2.855  '  10  "  ft. 

d  2 
^  =2. 125  btu  /(ft.  ^  sec. ) 

a  P 


A  bubble  initiation  diagram  of  these  points  is  shown  in  Fig.  27. 

A  comparison  of  experimental  results  with  theoretical  ones 
is  shown  in  Fig.  28 . 

5.  Description  of  Apparatus  and  Method  of  Experimentation 
a.  Experimental  Set-Up 

The  experimental  set-up  is  shown  in  Fig.  29.  The  heating 
surface  was  made  by  electroplat/nja  layer  of  16  k  gold  of 
0.005  inch  thickness  on  the  top  surface  of  a  thin  flanged 
cylindrical  copper  block.  The  reason  for  gold  plating  was  to 
minimize  the  effects  of  oxidation  so  that  the  surface  conditions 
will  remain  the  same  from  the  beginning  to  the  end  of  each  test. 
At  the  bottom  of  copper  block,  seven  120  watt  chromelux 
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electrical  heaters  were  imbedded  in  holes  in  the  copper 
block.  The  heat  generated  by  these  heaters  was  transferred 
to  the  top  surface  by  pure  conduction.  The  reduction  of 
cross  section  of  copper  block  underneath  the  heating  surface 
was  for  the  purpose  of  intensifying  the  heat  flux  at  the  heating 
surface.  A  thin  flange  surrounds  the  heater  to  eliminate  un¬ 
desired  bubble  nucleation  which  might  occur  at  a  boundary. 
This  flange  was  very  thin  so  that  the  temperature  near  the 
edge  of  the  heating  surface  was  low  enough  to  prevent  bubble 
initiation.  A  piece  of  Teflon  heat  insulator  was  inserted  be¬ 
tween  the  lower  face  of  this  thin  flange  and  pool  base.  A 
detailed  drawing  of  the  heating  surface  and  the  shank  part  of 
copper  block  is  shown  in  Fig.  30. 

A  thermo-bottle  filled  with  ice  was  used  for  the  cold  junction 
of  the  thermocouples  which  were  connected  with  a  potentio¬ 
meter  through  a  six -way  switch.  A  drain  hole  valve  was  also 
attached  to  the  bottom  of  the  test  section.  In  order  to  predict 
the  surface  temperature,  three  thermocouples  T^,  and 
T  were  inserted  in  the  holes  on  the  shank  part  of  copper 

u 

block,  a  three  point  interpolation  formula  was  used  to  deter¬ 
mine  the  wall  temperature  T  .  These  thermal  couple  holes 
were  1/16  inch  in  diameter,  19/32  inch  in  depth  and  were 
spaced  1/4  inch  apart.  All  dimensions  were  measured  from 
the  heating  surface.  The  bottoms  of  these  three  holes  were 
at  the  center  line  of  the  shank.  In  the  fluid,  another  thermo¬ 
couple,  T^,  was  used  to  measure  the  temperature  of  main 
body  of  fluid,  T^  .  It  was  located  one  inch  above  the  heating 
surface.  All  thermocouples  were  made  of  No.  30  Chrome - 
Alumel  wire.  In  order  to  avoid  excessive  corrosion,  the 
thermocouple  T^  was  shielded  in  a  1/16  inch  stainless  steel 
tube  with  Teflon  seal  at  the  outer  end. 

The  fluid  was  contained  in  a  3  inch  diameter  and  20  inch 
length,  specially  heat-treated,  high  strength  glass  tube. 
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Observations  and  photographs  could  be  made  through  a  so- 
called  "fluid  crystal".  This  was  a  glass  box  filled  with  the 
same  fluid  as  that  in  testing  section  and  so  placed  as  to 
eliminate  distortion  due  to  curvature.  The  front  wall  was 
flat,  the  rear  wall  was  made  of  a  segment  of  circular  tube 
with  a  radius  of  curvature  just  equal  to  the  outside  radius 
of  the  testing  tube,  such  that  the  distortion  of  bubble  shape 
due  to  light  refraction  of  tube  was  eliminated.  With  this 
device,  an  accurate  measurement  of  bubble  dimension  could 
be  obtained  from  high  speed  photography. 

A  helically  wound  copper  tube  in  the  upper  part  of  the  testing 

tube  was  used  as  a  condenser.  The  saturation  temperature 

of  the  fluid  Tga^.  was  controlled  by  varying  the  system  pressure 

from  1  atmosphere  to  1/4  atmosphere  through  an  asperator 

vacuum  pump.  The  temperature  of  the  main  body  of  fluid  T^ 

was  controlled  by  varying  the  flow  rate  of  the  cooling  water 

through  a  cock.  The  wall  temperature  T  was  controlled  by 

w 

varying  the  electrical  power  of  the  heaters  through  a  variac. 
b.  Surface  Preparation 

Boiling  data  are  difficult  to  reproduce  due  to  changes  in  the 
surface  conditions.  There  are  two  ways  in  which  these 
changes  appear;  namely,  changes  due  to  contamination  and 
cavity  reactivation.  Contamination  can  be  eliminated  by  pro¬ 
per  choice  of  the  metal  for  the  heating  surface,  reactivation 
of  a  nucleate  cavity  can  not  be  eliminated  by  the  following 
method. 

The  16  k  gold  plated  surface  was  first  finished  by  200 
grit  emery  paper  which  was  continuously  wetted  by  a 
a  water  jet.  The  direction  of  stroke  was  kept  constant. 

The  surface  was  finished  by  stroking  in  one  direction 
till  all  scratches  were  eliminated  then  rotating  90°  to 
eliminate  all  the  scratches  in  the  other  direction.  The 
whole  piece  was  then  washed  in  a  water  jet.  Following 
exactly  the  same  procedure,  the  surface  was  finished  by 
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400  grit  and  600  grit  emery  paper.  The  surface  was 
then  cleaned  by  hot  water  jet,  alcohol  jet  and  hot  air 
jet  and  was  then  put  on  the  No.  4  diamond  compound 
wheel.  The  diamond  compound  should  be  put  on  the 
center  area  of  grinding  wheel  and  diluted  by  kerosene 
before  starting  grinding  operation.  The  piece  was 
held  gently  near  the  edge  area  of  wheel,  kerosene  was 
injected  on  the  wheel  clothe  occasionally.  Operation 
was  continued  until  the  scratches  due  to  600  grit  emery 
paper  disappeared  completely.  Then  the  piece  was 
taken  off  from  No.  4  diamond  compound  wheel,  the  hot 
water  jet,  the  alcohol  jet  and  the  hot  air  jet  were  then 
each  put  on  the  surface.  After  the  washing  process, 
the  piece  was  then  put  on  the  No.  6  diamond  compound 
wheel  and  then  No.  8  wheel  using  the  same  sequence  of 
operations  as  on  the  No.  4  wheel. 

An  unclean  piece  will  leave  some  dust  particles  on 
the  wheel  which  sometimes  make  some  unremovable 
scratches  on  the  surface.  To  make  a  good  surface, 
one  needs  usually  more  than  10  hours.  Patience  and 
cleanliness  are  the  two  most  important  characteristics 
of  a  surface  worker.  Scratches  due  to  the  grinding  com¬ 
pound  can  be  removed  only  by  its  next  number  grinding 
compound  as  recommended  here.  No.  4  diamond  com¬ 
pound  scratches  can  not  be  removed  by  No.  8  diamond 
compound  wheel  in  a  reasonable  length  of  time  without 
introducing  No.  6  diamond  compound  wheel. 

At  the  last  few  minutes  of  grinding  process  on  the  No.  8 
diamond  compound  wheel,  the  kerosene  jet  was  applied 
all  over  the  center  area  of  the  wheel,  such  that  the  dia¬ 
mond  compound  was  washed  to  a  very  dilute  condition, 
the  piece  was  then  put  near  the  center  part  of  the 
wheel  where  the  rubbing  speed  is  lower,  then  a  heavier 
pressure  was  applied.  After  one  to  two  minutes,  the 
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surface  would  become  shining,  mirror-like  ^smooth. 

It  was  washed  by  hot  water  jet,  alcohol  jet  and  hot 
air  jet,  it  was  then  introduced  in  the  pool  of  a  ultra¬ 
sonic  cleaner  for  2  minutes.  This  process  would  help 
to  wash  out  small  diamond  dust  particles  and  bits  of 
metal  which  were  trapped  in  the  cavities  on  the  surface. 
Then  the  surface  was  washed  again  by  alcohol  and  Methyl 
ether  jet.  The  surface  at  this  stage  was  assumed  to  be 
the  surface  required. 

After  each  test,  the  surface  was  renewed  by  going 
through  all  the  steps  immediately  after  No.  6  diamond 
compound  wheel.  It  needed  only  20  minutes  to  finish 
the  job. 

In  order  to  keep  surface  condition  unchanged,  every 
element  which  is  in  the  boiling  system  should  be  cleaned 
by  washing  soap,  hot  water  jet  and  distilled  water  jet 
before  each  test. 

Method  of  Experimentation 

After  making  a  new  surface  and  washing  all  the  parts,  they 
were  assembled,  distilled  water  was  introduced  into  the  top 
of  the  test  section.  Two  hours  of  vigorous  boiling  with  a 
moderate  heat  flux  was  maintained  for  degassing  purposes, 
then  the  heat  flux  was  reduced  until  there  were  no  active 
cavities  on  the  surface,  then  the  heat  flux  was  increased 
gradually  until  the  first  active  cavity  appeared  on  the  surface. 
This  was  the  starting  point  of  each  test.  A  steady  state  con¬ 
dition  was  assumed  to  be  reached  two  hours  after  the  heat 
flux  was  changed. 

During  each  run  the  following  measurements  were  made. 
Power,  fluid  temperature,  heater-unit  temperatures, 
system  pressure,  number  of  active  centers,  and  num¬ 
ber  of  new  sites  arising  from  the  change  in  heat  flux. 
Technically  the  later  are  called  the  initiated  cavities 
which  generate  bubbles  with  very  low  frequencies  such 


that  the  contribution  to  the  heat  transfer  is  negligible  . 

The  heat  transfer  to  the  fluid  through  the  heating  sur¬ 
face  was  determined  by  the  simple  conduction  formula 
knowing  the  temperature  gradient  in  the  shank  of  copper 
block.  The  wall  temperature  Tw  within  a  circular  area 
of  ]  3/16"  diameter  on  the  center  part  of  the  heating 
surface  was  assumed  to  be  uniform  . 

d.  Photographic  Technique 

High  speed  photographs  were  taken  with  a  Wollensack  camera. 

A  Kodak  Tri-X  negative  100'  film  for  high  speed  photography 
was  used.  About  2400  frames  per  second  was  taken  which 
necessitated  a  reduction  of  voltage  supplied  to  the  Wollensack 
camera  motor  to  about  70  volts  through  a  variac.  A  500  watt 
illuminating  lamp  was  installed  at  the  rear  of  the  test  section 
at  about  6  inches  away  from  the  tube  center,  so  that  the  heating 
surface  looked  shining  bright.  The  focus  of  the  camera  was 
very  carefully  adjusted  such  that  no  relative  motion  between 
the  circle  on  the  focusing  lens  and  the  bubble  to  be  photographed 
was  observed.  Each  two  marks  of  time  on  the  film  represent 
1/60  sec. 

The  camera  was  placed  as  close  as  possible  to  the  test  sec¬ 
tion  without  losing  the  sharp  focus  required.  A  reference 
wire  of  0.  040  inch  in  diameter  was  placed  besides  the  bubble 
which  was  to  be  photographed.  The  bubble  diameter  measure¬ 
ments  were  made  by  projection  on  a  microfilm  projector.  A 
geometric  mean  value  of  bubble  diameters  in  three  principal 
axis  directions  was  considered  as  the  bubble  diameter  for 
volume  calculation. 

e.  Temperature  Calibration,  Wall  Temperature  Prediction  and 

Heat  Flux  Determination _ 

i.  Temperature  Calibration 

For  this  special  kind  of  chrome -alumel  thermocouple, 
the  following  data  were  recorded. 
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Reading  at  boiling  point  of  water  =  5.  209  milli -volts 

Reading  at  melting  point  of  tin  =  12.  410  milli -volts 

Reading  at  freezing  point  of  water  =  0.  000  milli -volts 
Atmospheric  pressure  P  =  77.66  cm./-/j 
Boiling  point  of  water  at  P  =  77.  66  Cm  ^ 
is  from  T  =  100.  000  +  0.  03686  (P  -  76.  00)  - 
0.0000220  (P-  76.  00)2  =  100.061°C 

Melting  point  of  tin  at  atmospheric  pressure  =  231.89°C 

Freezing  point  of  water  at  atmospheric  pressure  = 

0.000°C 

From  Reference  (14),  a  three  point  interpolation  formula 
gives 

T  =  0.07267/  VIZ17.S4-0  -  Y)  °C  (89) 

Where  V  is  the  reading  of  thermocouple  from  potentio- 
meter  (  Yn.V.  ) 

T  is  the  corresponding  temperature  C°  . 

ii,  Wall  Temperature  Prediction 

Referring  to  the  sketch  of  heat  surface  and  heat  conduc¬ 
tor  shank  which  is  shown  in  Fig.  30,  the  following 
dimensions  were  obtained  by  an  accurate  measurement. 

S,  ■ —  o.  2  A-  2  " 

S,  =  0.23  V' 

=S  C  z  3  o" 

The  location  of  each  thermocouple  and  the  heating  surface 
can  be  described  by  coordinate  X's,  say 

-X'  —  O,  ^  =  S, 

=c,=  +  J  =0476''  Xw  =  5,  +  S.  +  S3  =  0.706" 

If  the  corresponding  temperature  sat  X^,  X 2,  X^,  and 

X  are  T1,  Tn,  T0,  and  T  ,  then  with  help  of  Ref.  (14), 
w  1  2  3  w 

the  wall  temperature  can  be  extrapolated  by  Lagrangian 
method  as 
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-r-  _  y-  -j-  (yW-X,y>  (Xty-Vj^)  7"  (Xw~*mJ 

/w“*  '*  ty-x*;  f  x,-x3>  4  Cx.-x^  cxa-xp  J  6Xj~*J(**-X*> 

substituting  the  value  of  X^,  X^,  X^,  and  into  above 
equation  yields 

Jw  =  .9‘7*Tx+WT9  (90) 

iii,  Heat  Flux  Determination 

From  conduction  equation,  one  has 


Q  =  As  kc  r—  ■=  As  K  -~7;' T*  - 

s  c  *X  X3 

A  =  cross  section  area  of  shank  of  conductor  = 

S  2 
0.00769  12  ft. 

Kc  =  conductivity  of  copper  =  219  btu/(°F  hr.  ft.)  = 

0.  060833  btu/tft.*  °F-sec.) 

X3  =  0.476"  =  0.03967  ft. 

-  Tg  in  degree  Fahrenheit 

Q  =0.0118  (T1  -Tg)  btu/sec.  (91) 


6.  Discussion  and  Conclusion 
a.  Discussion 

In  the  preceding  sections  it  has  been  shown  possible  to  con¬ 
sider  the  individual  processes  of  bubble  initiation,  growth, 
and  departure  and  with  nothing  other  than  geometric  ideali¬ 
zations  and  fluid  and  surface  properties,  compute  a  heat  flux 
versus  wall  temperature  curve.  The  computed  and  measured 
heat  flux  curves  compare  satisfactorily.  In  making  this  com¬ 
parison  however,  an  extraordinary  amount  of  information  was 
needed.  In  practical  terms,  quantities  like  surface  nucleation, 
properties  and  bulk  temperatures  are  just  not  known  with  suf¬ 
ficient  precision  to  make  a  boiling  curve  prediction  possible. 

In  addition,  only  the  isolated  bubble  portion  of  the  pool  boiling 
curve  has  been  studied.  If  the  forced  convection  or  the  close 
packed  regions  are  of  interest,  then  other  geometrical  ideali- 
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zations  of  the  fluid  mechanics  are  needed.  This  all  raises 
the  question  of  where  do  we  go  from  here  ? 

The  close  correlation  between  theory  and  experiment  and 
the  fact  that  no  arbitrary  constants  have  been  used  show  that 
no  important  physics  has  been  forgotten.  What  further  should 
be  done  in  nucleate  boiling  ?  First,  the  physics. 

i.  Bubble  Initiation 

More  experience  is  needed  in  the  experimental  control 
of  surface  conditions.  Contact  angle,  drift  and  surface 
nucleation  properties  need  more  study  in  order  that 
meaningful  experiments  can  be  desired. 

ii.  Bubble  Growth 

This  appears  to  be  well  understood 

iii.  Bubble  Departure 

The  zero  gravity  departure  prediction  should  be  tested. 
Forced  convection  bubble  departure  should  be  studied 
as  only  a  few  odd  measurements  now  exist. 

Engineering  work  -  better  heat  flux  temperature  difference 
correlations  are  now  possible  based  on  our  cleaner  under¬ 
standing  of  the  basis  processes.  Work  with  industrial 
rather  than  laboratory  type  data  is  probably  most  desirable. 
Conclusions 

i.  The  nucleate  pool  boiling  curve  in  the  isolated  bubble  re¬ 
gion  can  be  predicted  from  a  knowledge  of  fluid  properties 
and  surface  conditions.  Resort  need  not  be  made  to  any 
physically  unmotivated  quantities. 

ii.  Dynamic  effects  on  bubble  departure  size  manifest  them¬ 
selves  primarily  through  contact  angle  variations  with 
the  usual  Fritz  formula  still  holding. 

iii,  A  formula  including  dynamic  effects  in  the  liquid  has  been 
developed  that  would  predict  bubble  departure  at  zero 
gravity  for  certain  fluid  properties  and  temperature 
distributions . 
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iv,  The  waiting  period  between  bubbles  is  shown  to  be  a 
known  function  of  cavity  size  and  liquid  and  surface 
temperature. 

v.  Using  measured  delay  times,  bubble  growth  rate  can  be 
predicted  with  good  precision. 

vi.  Using  measured  contact  angles  bubble  departure  size 
can  be  predicted. 

vii.  Contact  angle  has  been  found  to  be  a  function  of  velocity 
across  the  surface.  This  in  turn  has  been  correlated 
with  viscous  effects  in  terms  of  a  ratio  of  Webber  to 
Reynolds  numbers. 
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BUBBLE  GROWTH  TABLE  1 


The  Experimental  Data  for  Bubble  Number  1 

t  =  0.0245  sec.  Camera  speed  = 

w  1140  frames/sec. 

t ,  =0.  0166  sec. 

d 

R,  =  3.974  •  10_3  ft. 
d 


No. 

of 

Frame 

t 

Millisec 

Bubble  Diameter  on 
Microfilm  Projector 
(m .  m )  Scale  =  8.87*  1 

R 

Millift 

1 

0 

2.  32 

0.429 

2 

0.877 

9.48 

1. 754 

3 

1.  654 

13.  72 

2.  538 

4 

2.  631 

16.04 

2.  967 

5 

3.  508 

17.80 

3.  293 

6 

4.  385 

18.40 

3.404 

7 

5.  262 

18.95 

3.506 

8 

6.  139 

19.57 

3.620 

9 

7.016 

19.86 

3.674 

10 

7.893 

20.00 

3.700 

11 

8.770 

20.36 

3.767 

12 

9.647 

20.  72 

3.833 

13 

10. 524 

21.  10 

3.904 

14 

11.401 

21.  65 

4.005 

15 

12. 278 

21.90 

4.052 

16 

13.155 

21.90 

4.052 

17 

14.032 

21.80 

4.033 

18 

14.909 

21.69 

4.013 

19 

15.786 

21.48 

3.974(Rd) 
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BUBBLE  GROWTH  TABLE  2 


The  Experimental  Data  for  Bubble  Number  2 

t  =  0.0437  sec.  Camera  speed  = 

w  1260  frames/sec. 

t,  =0.0167  sec.  ' 

d 

R  ,  =  5.  328  *  10”3  ft. 
d 


No. 

of 

Frame 

t 

Millisec 

Bubble  Diameter  on 
Microfilm  Projector 
(m.  m  )  Scale  =  8 . 87  *  1 

R 

Millift. 

1 

0.  793 

9.26 

1.713 

2 

1. 586 

14.22 

2.631 

3 

2.379 

17.53 

3.243 

4 

3.  172 

19.  76 

3.656 

5 

3.965 

21.65 

4.005 

6 

4.758 

22.80 

4.218 

7 

5.551 

23.87 

4.416 

8 

6.344 

24.  67 

4.564 

9 

7.  137 

25.  20 

4.662 

10 

7.930 

25.48 

4.  714 

11 

8.723 

25.  50 

4.718 

12 

9.516 

25.80 

4.773 

13 

10.309 

26.  10 

4.829 

14 

11. 102 

26.  74 

4.947 

15 

11.895 

27.  26 

5.  043 

16 

12. 688 

27.40 

5.069 

17 

13.481 

27.  70 

5. 125 

18 

14. 274 

27.90 

5. 162 

19 

15.067 

28.53 

5.278 

20 

15.860 

28.87 

5.341 

21 

16.653 

28.80 

5.328(Rd) 
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BUBBLE  GROWTH  TABLE  3 


The  Experimental  Data  for  Bubble  Number  3 

t  =  0.0275  sec.  Camera  speed  = 

1380  frames/sec. 

t^  =  0. 0145  sec. 

Rd  =  0.395  '  10"3  ft. 


No.  ,  Bubble  Diameter  on 


of 

Frame 

Millisec 

Microfilm  Projector 
(m.  m  )  Scale  =  8. 87  :  1 

n 

Millift. 

1 

0.  725 

9.52 

1.711 

2 

1.450 

14.  10 

2.609 

3 

2.  175 

16.88 

3.  123 

4 

2.  900 

18.45 

3.413 

5 

3.  625 

19.55 

3.617 

6 

4.350 

20.  60 

3.811 

7 

5.075 

21.  73 

4.020 

8 

5.800 

22.08 

4.085 

9 

6.525 

22.53 

4.  168 

10 

7.  250 

22.  75 

4.  209 

11 

7.975 

23.  68 

4.381 

12 

8.  700 

24.  00 

4.440 

13 

9.425 

24.  10 

4.459 

14 

10. 150 

24.  33 

4.  501 

15 

10.875 

24.  15 

4.468 

16 

11. 600 

24.  74 

4.577 

17 

12.325 

25.  52 

4.  721 

18 

13. 050 

25.80 

4.  773 

19 

13. 775 

25.88 

4.788 

20 

14.500 

25.90 

4.792(Rd) 
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TABLE  4 


HISTORY  OF  BUBBLE  GENERATIONS 


T  =  229.98°F,  T  =  212°  F,  T  =  205.02°F 
w  sat  co 

Distilled  water  on  gold  surface  ground  by  No.  8  diamond  compound 


Bubble  Camera  Speed  t  (sec.)  t,  (sec.)  R ,  (Millift. ) 
No.  (frames /sec.)  w  d  d 


1 

1140 

0.0254 

0.0167 

3.974 

2 

1260 

0.0436 

0.0167 

5.328 

3 

1380 

0.0275 

0.0145 

4.  792 

4 

1500 

0.0466 

0.0167 

3.  534 

5 

1650 

0.0735 

0.0261 

3.  691 

6 

1800 

0.0594 

0.0172 

4.  224 

7 

1920 

0.0490 

0.0151 

4.  188 

8 

2040 

0.0633 

0.0162 

4.  658 

9 

2130 

0.0319 

0.0155 

3.  931 

10 

2190 

0.0337 

0.0160 

5.  125 

11 

2280 

0.0672 

0.0149 

3.321 

12 

2370 

0.0785 

0.0167 

3.448 

13 

2520 

0. 1640 

0.0143 

3.633 

14 

2700 

0. 1250 

/  0.0512 
*Three  in 
Tandem) 

4.  201 

15 

2770 

0.0436 

0.0143 

3.  566 

16 

28  50 

0. 0393 

0.0161 

4.  782 

17 

2910 

0.0216 

0. 0158 

5.  367 

18 

2910 

0. 0450 

0.0139 

3.  374 

19 

2940 

0.0756 

,  0.0296 
'Two  in 
Tdndem) 

3.  571 

20 

2940 

0.0252 

0.0163 

4.967 

21 

2940 

0.0354 

0.0139 

3.883 

22 

2940 

0.0490 

0.0129 

4.  183 

Observation  from  above  table  shows  that  the  waiting  period  t  changes 

from  17(t  )  ,  to  130(t  ) 

w'min.  w  mm. 
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meri  jqs  in  -them 


Fig.  1  Analogues  Between  Bubble  Initiation  &  Egg  Incubation 


Fig.  2  Temperature  Distribution  in  a  Semi-Infinite  Conductor 


Fig.  3  Transient  Thermal  Layer 
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X 


Fig.  4  Temperatures  of  Fluid  and  Bubble  Near  a 
Heating  Surface 


X 


Fig.  5  Initiation  of  Bubble  Growth  From  Different  Cavities 
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Fig.  11  Bubble  Growth  Plot 
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Bubble-  Number  /.  C.ame.ra,  Spa  Got  =  /  /40  -frames/ sec. 

Wai -f-mj  period  =  2.J  -frames  —  °PZ^^ 

tDepat  +ure-  Bgt! od  —  /  ?  freme S  =  otoj 66 sec 

t^GpQr-hjr-e-  Radius  wmm‘  3.  J!  Jo*  ft 


Fig.  12  History  of  ^Growth  of  Bubble  Number  3 
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Bubble  Mumbe-r 


Come.ro,  Spe&d  =  iZbo  frames  /soc 

Woi’/’tnj  period  =  £•£  frames  =  0.0+37  Sec. 

Vepar-bureu  period  ==  2./  frames  =<*0/67  Sec. 

t>*f>ar±ure.  Rod/uS  —  £.33  /o3  P 

Fig.  13  History  of  Growth  of  Bubble  Number  2 
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Bubble.  A/ umber  S. 


Camera.  Speed  =  /38o  promos  /sec 
Wa:iin<j  period  =  3g  pta„eS=  0.027*  Sec 
PeparbUr^.  Period  2.o  Frames  Sec. 

peparFure,  Radius  ==■  4.7?  /o3  Ft 

Fig.  14  History  of  Growth  of  Bubble  Number  1 
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Fig.  10  Normalized  Bubble  Growth  Curves 


Fig.  16  Dynamic  Load  on  a  Growing  Bubble 
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Fig.  17  Staniszewski's  Departure  Criterion 
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Fig.  18  Bubble  Generating  Cycles 


0,2$ 
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/  lo  /o 2  /o*  fo* 

7*'T-  (°p> 


Fig.  20  Boiling  Curve 


Fig. 21  Physical  Model  of  Bulk  Convection  Mechanism 
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Fig.  24  Nucleate  Cells  at  Close  Packed  Condition 


Fig.  25  Scavenging  Effect  of  a  Departing  Bubble 
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Wall  diS-htnCG.  ( /cT^Pf  ) 


Fluid,  to  “ 

Surface...  A//cke.l  ,  4/o  Tohsh^d 

~Too  88  TSa+ 

P  =•  /  a*//* 


Fig.  27  Bubble  Initiation  Diagram  of  Corty  and  Foust's  Data 
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Fig.  29  Pool  Boiling  Apparatus 
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Fig.  30  Detail  of  Heating  Surface  and  Heat  Conductor  Shank 
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APPENDIX 


THE  EFFECT  OF  VELOCITY  ON  CONTACT  ANGLE 


It  is  evident  from  Figures  12,  13,  and  14,  that  the  contact 
angle  is  smallest  when  the  contact  ring  is  growing  most  rapidly 
and  largest  when  the  contact  ring  is  contracting  most  rapidly.  The 
data  reported  here  are  too  skimpy  to  get  much  of  a  measure  of  the 
importance  of  this  effect,  so  the  purpose  of  this  appendix  is  to  re¬ 
late  the  deviations  in  contact  angle  from  the  static  condition  to  the 
growth  velocity  of  the  bubble  which  is  a  known  quantity.  The  ex¬ 
tensive  data  of  Staniszewski  will  then  be  used  to  determine  the  best 
value  for  the  constant  in  the  relationship  between  these  two  quantities. 

The  relationship  between  the  apparent  contact  angle  and  con¬ 
tact  ring  radius  is  shown  in  Figure  31  for  the  three  bubbles  illustrated 
in  Figures  12,  13,  and  14.  For  these  three  bubbles,  the  data  of 
Figure  31  for  the  rate  of  contact  ring  radius  change  can  be  obtained 
at  the  point  of  departure  along  with  the  associated  contact  angle. 

These  points  are  plotted  on  Figure  32  as  three  circled  crosses.  At 
the  same  time,  the  growth  rate  of  the  bubbles  at  departure  can  be 
plotted  versus  the  contact  angle.  These  points  are  the  circled  X’s. 
There  is  obviously  a  similar  relationship  between  these  two  velocities 
and  the  contact  angle  deviations.  As  the  growth  velocity  is  a  known 
quantity  in  this  work,  it  is  this  quantity  which  has  been  used  to  cor¬ 
relate  rather  than  contact  ring  velocity. 

To  repeat,  the  physical  velocity  of  significance  in  determing 
contact  angle  deviations  is  the  contact  ring  velocity.  This,  in  general, 
is  not  known  but  bears  a  casual  relationship  to  the  growth  velocity,  so 
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it  is  this  quantity  used  to  correlate  the  contact  angle  deviations. 
The  constant  in  the  equation  relating  the  static  contact  angle  to  the 
dynamic  contact  angel  equation  (64)  is  determined  from  the  exten¬ 
sive  data  of  Staniszewski.  Figure  32  gives  the  calculated  and 
observed  values  of  all  these  quantities. 
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NOMENCLATURE 


Dimensions  in  H,  M,  L,  T,  9;  The  Heat  Energy,  Mass,  Length,  Time, 
and  Temperature. 


A 

D 

K 

K 


L 

Mi 

N 

Na 

N. 


1 

P 


QR 

Qp 


Area  of  heating  surface 

Surface  characteristic  length  for  natural  convection 
Thermal  conductivity  of  fluid 
Thermal  conductivity  of  copper 


[L] 

[Hr'Ce" } 
[ht'l-9"] 


Latent  heat  of  evaporation  of  fluid  [H  M-J 

Inertia  mass  of  fluid  surrounding  a  bubble  [  M] 

Total  number  of  nucleate  center  son  heating  surface  [  O  ] 

Total  number  of  active  nucleate  centers  on  heating  surface  [  O  ] 


Total  number  of  initiated  nucleate  centers  on  heating  surface  [  O  ] 


Pressure  in  the  fluid  outside  the  bubble 
Heat  flux  received  by  heating  surface 

Heat  flux  predicted  by  theory 


[ML"TaJ 

[HT"] 

IHT-J 


R  Radius  of  bubble 

Rc  Radius  of  cavity 

R^  Departure  radius  of  bubble 

R  Meridian  curvature  radius  of  bubble  at  the  base  circle 
S  Bubble  surface 

T  Temperature 

T^  Temperature  of  vapor  in  the  bubble 

Tgat  Saturation  temperature  of  fluid  at  system  pressure 
Tjv  Wall  Temperature 
T^  Temperature  of  main  body  of  fluid 

u  Rising  velocity  of  a  solid  sphere 
v  Reading  of  thermocouple 


[L] 

[L] 

[L] 

[L] 

[L2] 

[6] 

[9] 

[9] 

[9] 

[0] 

[LT"1] 
[Milli  volts] 
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a 

b 

c 

f 

g 

h 

A 


K 

it 

n 

n« 

71- 

P 

} 

A. 

t 

t, 

tk 

tw 

X 

1 

A 

3> 

* 

s 

n 

e 

A 

s 


Radius  of  a  solid  sphere 

Radius  of  curvature  of  bubble  at  its  top 

Specific  heat  of  fluid 

Frequency  of  bubble  generation 

Gravity  acceleration 

Distance  from  the  center  of  a  solid  sphere  to  the  solid 
plane  boundary 

Coefficient  of  heat  transfer  from  wall  to  the  fluid 

Coefficient  of  heat  transfer  from  wall  to  Vapor 

Thermal  difusivity  of  fluid 

Number  of  nucleate  centers  per  unit  area 

Number  of  active  nucleater  center  per  unit  area 

Number  of  initiative  nucleater  center  per  unit  area 

Pressure  inside  the  bubble 

Heat  flux  density 

Normalized  bubble  radius 

Time 

Departure  period 
Unbinding  period 
Waiting  period 

Distance  from  the  bubble  surface  to  the  axis  of  revolution 
of  bubble;  Wall  distance 

Distance  from  the  plane  tangent  to  the  bubble  at  top  paint  of 
bubble  to  a  point  on  the  bubble  surface 

Constant  in  dynamic  contact  angle 

Velocity  potential 

Stream  line  function 

Volumetric  thermal  expansion  coefficient  of  fluid 
Thermal  layer  thickness 

Distance  from  the  bubble  center  to  base  plane  angle 

T'Tsat  ;  Anl'* 

Coefficient  of  Viscosity 

Kinematic  Viscosity 

Density  of  fluid 

Density  of  Vapor 


[L] 

[L] 

(wM-'e'1  J 

[  T-1] 

U  T2] 

ID 

[htVs"  I 
U*t"  I 

rm 

[t-n 

il-‘] 

[ML-'T2] 

[Hr'r1] 

[O] 

[0] 

[0] 

[0] 

[0] 

[L] 

[L] 

[  O] 

U,t"J 

re"i 

[i-] 

[0];  fO] 

[M  ] 

IM  L~3  ] 

[M  Cs] 
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^  Surface  tension  of  fluid 
Normalized  time  variable 
Volume  of  bubble 

Angle  of  Contact  in  Static  Condition 

Dynamic  contact  angle  at  the  instant  of  bubble  departure 
Base  factor 
Curvature  factor 
Surface  factor 
Volume  factor 
Nusselt  number 
Rayleigh  number 


LM  f1] 

[O  ] 

[L3] 

[O] 

[  O] 
[O] 
[O] 
[O] 
[O] 
[O] 
[O] 
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Subscripts 

Be  Bulk  convection 
cp  Close  packed  condition 
^  Departure 
Me  Natural  convection 
Saturation 
Ub  Unbinding 
Vc.  Vapor  convection 
w  Wall;  Waiting 
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